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Preface 


The  mercury-chloride  laser  is  a potentially  outstanding 
high  energy  laser.  The  need  for  understanding  the  signifi- 
cant microscopic  processes  is  essential  if  a large  scale 
system  is  to  be  built. 

This  paper  is  an  attempt  to  develop  the  temporal 
evolution  of  the  mercury-chloride  discharge.  A complete 
numerical  solution  of  the  electron  Boltzmann  equation  is 
developed.  Cross-sections  for  higher  metastable  levels  of 
mercury  are  determined.  A numerical  solution  for  the 
change  of  each  species  is  incorporated  with  the  numerical 
Boltzmann  code  for  the  overall  temporal  evolution. 

Several  people  have  provided  me  with  assistance  and 
aid  throughout  this  study.  Without  the  exchange  of  ideas 
I have  had  with  these  people,  this  paper  would  lack  some 
key  elements.  In  particular,  I want  to  thank  Ed  Seward  of 
my  GEP  78-D  class  for  the  numerous  hours  of  assistance  on 
the  computer.  I want  to  thank  Hal  Hastings  and  Tom  Gist 
for  the  joint  effort  in  numerically  solving  the  Boltzmann 
equation.  A large  debt  of  gratitude  is  owed  to  A.  M.  Hunter 
of  the  AFIT  Physics  Department,  my  thesis  advisor,  for  the 
guidance  and  encouragement  when  things  looked  as  though 
they  would  fall  apart.  I want  to  thank  Mrs.  Sharon  Maruna 
for  doing  a fine  job  in  typing  the  final  copy  of  this 
thesis.  Finally,  a special  note  of  thanks  to  my  wife  Kathy, 
for  her  patience  and  understanding. 
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Abstract 

Computational  simulations  of  the  HgCl  discharge  are 
discussed.  Equations  that  model  the  laser  are  solved 
simultaneously  as  functions  of  time.  The  discharge  model 
includes  electron  losses  by  dissociative  attachment  and 
dissociative  recombination.  Electron  sources  are  direct 
electron  beam  ionization  and  ionization  from  the  excita- 
tion of  the  atomic  species.  The  excitation  and  ionization 
rates  are  obtained  from  solutions  to  the  electron  Boltzmann 
equation  which  includes  the  following  processes:  inelastic 
collisions  of  the  first  and  second  kind,  heating  by  the 
electric  field,  secondary  electron  production  by  the 
electron  beam,  and  electron-electron  collisions. 

The  HgCl  discharge  model  determines  two  conclusions. 
First,  the  significant  metastable  levels  of  Mercury  for 
determining  the  reaction  rate  constants  are  the  6^ p . and 
6 p^  levels.  Second,  that  the  discharge  is  highly  sensitive 
to  the  concentration  of  chlorine. 
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TEMPORAL  EVC„  UTION  OF  THE  MERCURY - 


CHLORIDE  DISCHARGE 
I . Introduction 

The  mercury-chloride  molecule  is  presently  being 
studied  as  a possible  source  for  an  efficient,  high  energy, 
visible  laser.  Lasing  action  in  HgCl  from  an  electron  beam 
controlled  discharge  was  reported  in  August  of  1977  by  a 
joint  Maxwell-SRI  team  of  Tang,  Hunter,  Oldenettel,  Hov/ton, 
Huestis,  Eckstrom,  Perry,  and  McCusker  (Ref  1).  In  addition, 
Parks  (Ref  2)  and  Eden  (Ref  3)  report  observing  lasing  in 

O 

the  mercury-chloride  laser  at  a wavelength  of  5576 A. 

As  the  mercury-chloride  laser  is  developed,  the  need 
grows  to  understand  better  the  microscopic  processes  which 
determine  the  population  of  the  upper  and  lower  laser  levels. 
Understanding  these  processes  will  enable  experimentalists 
to  design  a variable  large-scale  system.  The  best  method  to 
understand  the  system  is  to  determine  all  the  important  con- 
tributing reactions  theoretically  within  an  order-of-mag- 
nitude.  These  reactions,  or  their  corresponding  cross 
sections,  can  be  compared  with  experimental  data  to  deter- 
mine both  their  importance  and  exact  value.  In  this  paper, 
as  many  reactions  and  cross  sections  as  can  be  determined 
significant  are  tabulated.  This  information  is  then  in- 
corporated in  a numerical  program  to  determine  the  temporal 
evolution  of  the  mercury-chloride  discharge. 
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The  mercury-chloride  lasing  action  is  obtained  by 


using  an  electron-beam  controlled  discharge  to  pump  a gas 
mixture  of  mercury-chlorine-argon.  A microdensitometer 
trace  of  the  lasing  action  is  given  in  the  article  by 
Hunter,  et  al.  (Ref  1)  and  is  shown  in  Figure  1. 


Fig.  1.  Microdensitometer  trace  of  B X HgCl  laser 

transition:  shown  here  is  also  the  reference 

Hg  5461 -A0  line.  Pressure:  Hg/Cl?/Ar  - 
10/8/1500  torr  respectively.  (Ref  1). 

In  addition  to  the  obvious  characteristic  of  lasing 
in  the  visible  spectrum,  the  mercury-chloride  laser  has  two 
other  important  features.  The  potential  quantum  efficiency 
is  4 and  the  shorter  wavelength  means  more  energy  per 
photon  than  the  infrared  lasers.  The  quantum  efficiency  is 
determined  by  the  ratio  of  the  emitted  photon  energy  to  the 
energy  necessary  for  creating  the  upper  laser  level. 
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The  approach  of  this  paper  is  to  analyze  the  mercury- 
chloride  gas  evolution  in  a three  step  sequence.  First,  a 
solution  of  the  electron  Boltzmann  equation  is  constructed 
to  determine  the  distribution  of  electrons  in  energy  space. 
Second,  the  metastable  levels  of  mercury  are  determined 
along  with  their  appropriate  cross  sections.  These  cross 
sections  are  necessary  for  calculating  metasble  and  ioniza- 
tion pumping  rates  by  the  electron  Boltzmann  equation.  In 
the  third  step,  the  physics  for  determining  the  temporal 
evolution  of  the  mercury-chloride  discharge  is  developed. 

The  third  step  includes  a tabulation  of  all  reactions 
thought  to  apply  to  the  mercury-chloride  laser.  The  reac- 
tions are  an  extension  of  the  work  of  Summers  (Ref  4)  and 
include  the  following  types: 

1 . electron-beam  ionization 

2.  electron  loss  reactions 

3.  metastable  formation 

4.  secondary  ionization 

5*  excimer  formation 

6.  excimer  loss  mechanisms 

7.  three  body  reactions 

8.  Penning  ionization 

9.  lower  level  eximer  loss 

10.  metastable  transfers 

Any  assumptions  made  in  this  paper  are  included  with 
the  establishment  of  the  appropriate  theory.  The  validity 
of  classical  mechanics  and  quantum  mechanics  is  a general 
assumption  made  throughout  the  paper. 

The  following  notation  is  used  when  referring  to 
specific  species; 


Mercury-Chloride  (lower  laser  level) 

Argon 

Mercury 

Chlorine 

Argon  (First  Metastable) 

Argon  (Second  Metastable) 

Mercury  (First  Metastable) 

Mercury  (Second  Metastable) 

Mercury  (Third  Metastable) 

Mercury  (Fourth  Metastable) 

Argon  (ionized) 

Mercury  (ionized) 

Mercury-Chloride  (upper  laser  level) 


HgCl 

Ar 

Hg 

Cl2 


Ar* 

Ar** 


Hg*** 
Hg 
Ar+ 
Hg+  , 
HgCl* 


The  term  "pumping  rate"  refers  to  the  rate  at  which  a 
species  is  excited  into  a metastable  level  by  a secondary 
electron.  Reaction  rate  constants  are  those  rate  constants 
for  reactions  other  than  a pumping  rate.  Any  time  deriva- 
tive, dn,  is  abbreviated  by  the  more  common  form  n and  n.  is 
dt  1 

the  abbreviation  for  the  number  density  of  species  i. 
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II . Boltzmann  Equation 


Introduction 

In  this  chapter,  the  electron  energy  distribution  in 
a plasma  with  an  applied  steady  electric  field  is  investi- 
gated. A simple  analysis  to  the  electron  Boltzmann  equa- 
tion is  developed  and  then  expanded  to  include  necessary 
terms  for  an  accurate  numerical  solution.  The  complete  ex- 
planation of  the  numerical  solution  is  maintained  by  Hunter 
(Ref  21).  Results  of  the  numerical  solution  are  compared 
to  results  of  Rockwood  (Ref  5)  and  results  by  a numericax 
code  supplied  by  Stamm  (Ref  29). 

Boltzmann  Equation 

Analysis  of  the  electron  energy  distribution  is  found 
throughout  the  literature  (Ref  5»  7,  and  8).  The  following 
development  follows  the  work  of  Nielsen  (Ref  7). 

The  basic  function  governing  the  electron  energy  dis- 
tribution is  the  Boltzmann  equation  shown  in  equation  1 . 


+ v-  r f(v)+  5-  vyf(v)  = ff^-j  (1) 

' Collisions 

The  first  three  terms  of  equation  1 represent  re- 
spectively the  change  of  the  velocity  distributi.  n with 
respect  to  time,  the  diffusion  of  particles,  and  the  motion 
imparted  by  an  external  force. 
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For  a spatially  uniform  gas  with  a constant  electric 


field,  equation  1 reduces  to  equation  2. 


SfjJD  + 
at 


a 


* vf<*> 


3l~(  v) 
at 


Collisions 


(2) 


In  many  cases  of  physical  interest,  the  second  term 
on  the  left  hand  side  of  equation  2 is  not  expected  to  im- 
part much  net  motion.  Therefore,  the  electron  distribution 
is  considered  symmetric  with  a small  anisotropy  in  the 
direction  of  v.  This  prompts  the  expansion  of  f(v)  as 


f(v)  = fQ(v)  + f1(v)  cos  e (3) 

where  0 is  the  angle  between  the  applied  field  and  v. 

Expressing  v in  terms  of  spherical  coordinates  (with 
a in  the  direction  of  the  polar  axis),  substituting  equation 
3 back  into  equation  2,  multiplying  equation  2 by  (cos  0)0'1 
then  integrating  over  all  solid  angles  yields 


3f0  + a ( L_  *1  ) 1 r Q d O 

+ ^ 1 > Err  j 


at  3 y2  av 


5fi  a afo  3 p cos  e q d n 

at  dv  4tT  J 


where  Q represents  the  collision  term. 


(4) 


(5) 
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Expressing  the  right  hand  side  of  equation  5 as  an 
effective  collision  time  (t)  for  f ^ , the  equation  can  be 
rewritten  as  equation  6. 


3fj  a afQ  f 
3t  + 3v  " r 


(6) 


The  solution  to  this  differential  equation  is 


f 


1 


A e 


-vt 


v 3v 


(?) 


where  v equals  */t. 

Ignoring  the  initial  transient  time,  and  substituting 
equation  7 into  equation  4,  the  result  is  equation  8. 


af0  a?  1_  3_  ["  v2  ^£o  1 1 . Q d o 

3t  ‘ v y2  3v  [ 3 v J J 


(8) 


In  equation  8,  v is  assumed  to  be  independent  of  v. 
Converting  to  the  more  useful  energy  space,  by  way  of 
equation  9 and  10, 


f(€)d€  = 4rTfQ(  v)  v2dv, 

(9) 

v(€)  = (2e/m)l/2, 

(10) 

and  recalling  a m = eE,  the  final  result  is  equation  11. 
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W\ 


Q d fl 


(11) 


3f  3J 
d t 9 £ 


2£ 


1/2 


m 


372* 


where  J = ~- 


2 2 
2 e^E  € 


3 mv 


f L_  _ af  1 

L 2€  d€  J 


Relating  f(€)  to  fQ(€)  "by  way  of 

4tt  i/2  €1//2  fn(€) 


f(€)  = 


m 


w 


(12) 


and  expressing  Q as  an  inelastic  collisional  process,  the 
right  hand  side  of  equation  11  becomes 


ng  f(€-A€)  R (e-A€)  - npf(€)  R(€) 


-np  f(e)  R'(e)  + ng  f(e  + aO  R(e  + Ae) 


(13) 


where  R(0  = o(£)  v(€)  is  the  rate  of  producing  electrons 
with  velocity  v(€),  losing  an  energy  a € , and 

R'(€)  = o'(€)  v(€)  is  the  rate  of  superelastic 
collisions.  Also, 


•’(e)  = [ — — A-  j 0 (e  + a€) 


by  way  of  microreversibility.  The  result  is  identical 
to  Rockwood's  derivation  (Ref  5)  without  electron-electron 
collisions  and  ionization  processes. 
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Numerical  Solution 


Determining  as  many  collisional  processes  as  possible 
is  necessary  in  developing  an  accurate  solution  to  the  elec- 
tron Boltzmann  equation.  The  entire  solution  is  maintained 
with  Hunter  (Ref  21),  and  includes  elastic  collisions  (in- 
cluding electron-electron),  inelastic  collisions,  heating, 
ionization  by  secondary  electrons,  and  ionization  due  to  an 
electron  beam. 

Using  a notation  similar  to  Rockwood  (Ref  5)i  the 
entire  Boltzmann  equation  is  written  in  equation  14. 


an  (e ) _ 
at 


3 J_  _ a J r.  J, 

f + R + R + E s_(e)  (14) 

s s 


ae 


ae 


where  Jee  = v [n  ( |^  - e ) - KT€  ] 


v = 2mN  (26/m)1/2  E qo  n (€)/M 


s s 


= 


2 N e2  (E/N)2  6 r n_  an 

3m  ( v/N)  L ^ ^ J 


26  36 


y. 

N 


< ir  > J °B  <e> 


1/2 


* 

R 


N ^ 

E N°[  R . (6  + A6)  + R’  .(6-A€)n(€-A6) 

s,j  b SJ  N 0 

s 


-Rsj(€)N(€)  - Rg.(€)N(€)  ] 
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• oo 

R+  = C K1  (e  + Ae)n(e  + a€)  + 6(e)  f r 1(€)n(e)d€ 

s A€  s 

-Rs1(e)n(e)  ] Ns° 

2 Sg(€)  = rate  at  which  an  external  electron  beam  creates 
secondary  electrons. 

In  equation  14,  q represents  the  mole  fraction  of 

D 

species  s,  a_  the  momentum  transfer  cross  section  from 
s 

electrons  of  energy  € to  molecules  N , e and  m are  the 

o 

electron  charge  and  mass.  T represents  the  gas  temperature, 

and  M the  mass  of  species  s.  R represents  the  rate  constant 

for  all  inelastic  collisions  of  the  first  and  second  kind  of 

species  s to  level  j,  while  R+  represents  the  rate  constant 

for  ionization  of  all  species  N . In  the  term  R+,  6(0  re- 

s 

presents  the  creation  of  secondary  electrons  at  zero  energy. 
This  assumption  allows  for  a conservation  of  particles. 

The  results  of  the  next  section  were  determined  with 
the  numerical  code.  The  technique  of  developing  the  code  is 
shown  in  Appendix  C. 

Results 

The  drift  velocity  for  the  electron  distribution  is 
computed  and  compared  to  those  calculated  by  Rockwood.  The 
results  shown  in  Figure  2 are  from  equation  15. 
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The  bend  at  the  lower  E/N  value  is  not  found  due  to 
the  coarse  bin  width  used  in  the  numerical  solution,  and  the 
fact  that  most  of  the  energy  goes  into  elastic  heating 
which  is  not  included  in  the  numerical  solution  at  this  time. 
Elastic  collisions  were  included  at  a later  time  and  the 
results  matched  throughout  the  E/N  range. 

A similar  numerical  generated  electron  distribution 
function  for  an  Argon  pure  gas  is  supplied  by  Stamm  and  Bailey 
(Ref  29)  and  a comparison  is  shown  in  Figure  3- 
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f(0  (eV 
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As  can  "be  seen  from  Figures  2 and  3,  the  results  of 
the  numerical  solution  to  the  electron  Boltzmann  equation 
agree  rather  well  with  the  other  two  numerical  solutions. 


Pumping  Rate  Constants 


The  major  use  of  the  electron  distribution  function 
is  to  calculate  pumping  rates  and  drift  velocity. 

The' numerical  solution  of  the  electron  Boltzmann 


equation  uses  cross  sections  to  determine  the  pumping  rate 
constants  by  way  of  equation  1 6. 


oo 


oo 

S f(0  d€ 


0 


(16) 


In  equation  16,  R represents  the  pumping  rate  constant 

o 

(cm-ysec)  from  level  x to  level  y,  Q is  the  cross  section 
(cm  ) for  that  transition,  v(€)  is  the  velocity  of  the 
electron,  and  f(£)  is  the  electron  energy  distribution 
function . 

Before  the  electron  energy  distribution  can  be  used  to 
calculate  pumping  rate  constants,  the  cross  sections  for 
the  excitation  and  ionization  of  each  gas  species  must  be 
determined.  The  next  chapter  is  devoted  to  the  calculation 
of  appropriate  cross  sections  for  the  Hg/Cl2/Ar  mixture. 


III.  Cross  Sections 


Overview 

In  this  chapter,  the  cross  sections  necessary  for 
calculating  pumping  rates  from  the  electron  Boltzmann  equa- 
tion are  developed.  The  two  major  processes  discussed  are 
those  of  electron  collisional  excitation  and  electron  col- 
lisional  ionization. 

Of  the  three  gases  in  the  HgCl  laser,  Hg/Cl2/Ar, 

Hg  cross  sections  are  analyzed  in  detail.  Cross  sections 
for  the  Argon  reactions  have  been  previously  determined. 

These  values  are  listed  in  the  Argon  section  of  this  chapter. 
The  section  for  the  cross  section  by  dissociative  attach- 
ment of  Cl2  treats  the  process  as  a narrow  function  oc- 
curring around  zero  energy. 

Some  of  the  cross  sections  determined  in  this  chapter 
are  a result  of  direct  experimental  measurement,  while 
others  are  the  result  of  solutions  to  theoretical  and  exper- 
imental equations.  The  next  section  is  the  development  of 
the  Hg  cross  sections. 

Mercury  (Hg) 

Introduction . A metastable  level  is  by  definition 
(Ref  8)  a state  of  an  atom  which  cannot  transition  to  a 
lower  state  with  the  emission  of  radiation,  and  correspond- 
ingly, cannot  be  reached  from  a lower  state  by  absorption. 
Various  transition  rules  prohibit  these  "forbidden 
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transitions".  If  this  transition  does  occur  by  some  other 
method  than  photon  absorbtion,  i.e.,  collisions  of  the  first 
kind,  the  result  is  an  extremely  long  lifetime  for  the  ex- 
cited atom.  In  some  cases,  lifetimes  of  milliseconds  are 
not  uncommon.  Transitions  caused  by  collisions  of  the  first 
kind  or  collisions  of  the  second  kind  are  prohibited  only 
by  the  Pauli  principle.  Because  of  this  type  of  interaction, 
evaluation  of  the  entire  energy  level  structure  is  important. 

An  atom  may  be  collisionally  excited  into  a metastable 
state,  and,  before  being  radiatively  deexcited,  collisionally 
reexcite  into  a higher  energy  level.  The  lifetime  of  this 
higher  level  may  be  short,  causing  the  atom  to  deexcite  back 
into  the  metastable  state.  This  process  can  occur  numerous 
times  before  the  atom  changes  its  condition  by  some  other 


process.  The  atom  behaves  as  if  it  is  storing  energy.  If 
a secondary  electron  collides  with  this  atom,  while  it  is 
in  the  higher  state,  the  atom  may  be  collisionally  stripped 
(ionized)  of  its  electron.  In  essence,  an  atom  capable  of 
being  in  a stored  energy  condition,  very  often  behaves  as 
being  in  a metastable  state  of  the  higher  energy. 

A typical  energy  level  diagram  with  spectral  lines 
for  the  element  Hg  is  borrowed  from  Herzberg  (Ref  8)  and 
illustrated  in  Figure  4. 
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Table  I lists  the  probabilities  of  17  various  transi- 
tions for  the  element  Hg. 

Table  I 

Transition  Probabilities 


Transition 

A(  ,i . k)  ( sec ) 

Transition 

A(  .i,k)  (sec) 

73P2- 

73S1 

0.44 

73S1  — 

C\J 

Ph 

cn 

M3 

1.1 

73Pj— 

73si 

0.18 

73S1— 

63  Px 

1.2 

73P  — 

73S1 

0.23 

73S1  — 

63p0 

0.52 

73P1 — 

7^0 

0.035 

73Sx  — 

61pl 

0.002 

71P1  — 

7^0 

0.18 

7^0- 

61pl 

0.29 

71P1  — 

73S1 

0.15 

7^0- 

63P1 

0.13 

73P1 — 

6ls0 

0.13 

63d2- 

63P2 

0.54 

71P1  — 

6^0 

0.37 

6\~ 

63p1 

3.30 

63d2- 

6lpl 

0.16 

O 4 

Transition  Probabilities  (in  units  of  10  Sec’  ). 
(Ref  9)  Limits  of  error  omitted. 


For  ease  of  understanding,  the  radiative  lifetimes 
are  computed  by, 

T (j,k)  = A-1U,k)  (17) 

and  displayed  in  Figure  5-  The  units  of  A are  sec'^. 

The  lifetimes  of  the  b-^p^  ^ q level  are  measured  by 
Panevkin  (Ref  10)  to  be  4 X 10~3,  1 X 10~3,  and  6 X 10“3 
sec  respectively.  Most  of  the  transitions  shown  in  Figure 
5 are  fast,  occurring  between  1 X and  1 X 10-^  seconds. 
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Investigation  of  the  energy  level  diagram  with  the 

associated  lifetimes  indicates  the  possibility  of  trapped 

3 3 11 

radiation  between  the  states  7 6^Pg  1 O’  ^ Pj. 6 Sq. 

and  63D 6-^p ^ The  high  transition  probability  and  the 

long  lifetimes  of  the  lower  states  are  good  indicators. 

Rockwood  (Ref  5)  and  Tang  (Ref  11)  recognize  the  need 
to  consider  higher  excited  levels  of  Hg  other  than  the 

3 

Pg  j q level.  These  higher  excited  levels  may  act  as 
metastable  states,  capable  of  being  ionized  by  secondary 
electrons . 

In  order  to  arrive  at  more  accurate  results,  the 
temporal  evolution  of  these  excited  levels  needs  to  be 
tracked.  The  next  section  develops  the  cross  sections 
necessary  for  following  these  levels.  The  notation  and 
corresponding  energy  level  is  tabulated  in  the  following 
list. 


Notation 

State 

Level  (eV) 

* 

Hg 

6\ 

5 . 46 

** 

Hg 

6lpj 

6.70 

Hg 

73S1 

7.73 

Hg 

6\ 

8.8  5 

The  reason  for  tracking  the  63P2  level  instead  of  all 
three  6^P  levels  is  mainly  for  simplicity  and  ease  of  cal- 
culations. A quantitative  argument  for  tracking  only  the 
6^p2  level  is  made  in  the  results  section  of  this  chapter. 
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Theory 

Electron  Impact  Excitation 

A pumping  rate  constant  is  only  as  good  as  the  cross 
section  which  determines  that  rate  constant.  Cross  section 
curves  can  be  calculated  with  quantum  mechanical  methods; 
however,  the  calculations  are  extremely  complicated  and  re- 
quire approximations  that  make  the  degree  of  accuracy  almost 
impossible.  An  appreciation  of  this  difficulty  is  evident  in 
the  work  of  Yavorsky  (Ref  12).  He  calculates,  quantum 
mechanically,  the  cross  section  for  the  6^P^—  6^P2  Hg 
transition  by  a free  electron. 

Fortunately,  Bauer  and  Bartkly  (Ref  13)  show  that 
classical  theory  is  often  adequate  for  determining  colli- 
sional  cross  sections.  This  technique  works  best  when  the 
free  electron  is  below  27  eV  (e  /hv»l). 

A classical  differential  cross  section  is  defined  by, 

a = d Qn/d€  (18) 

In  equation  18,  6 is  the  energy  (eV)  transferred  from  the 
incident  electron  to  the  orbital  electron,  Qn(cm  /eV)  re- 
presents the  total  cross  section  for  transition  n.  An  in- 
tegration over  the  limits  of  energy  needed  to  reach  a level 
(Un+i)  yields  a total  cross  section  as  in  equation  19-  The 
upper  limit  of  integration  is  determined  by  the  energy 
necessary  to  reach  the  next  highest  energy  level. 


(19) 


«n  “ ^ °n‘€>  d f 


n 


With  the  approach  of  equation  19*  Bauer  and  Bartkly 
(Ref  13)  examine  several  theoretical  and  empirical  solutions 
for  the  differential  cross  section  g(0.  In  particular, 
they  cite  the  work  of  Thompson  (Ref  14),  Drawin  (Ref  15)*  a 
series  of  papers  by  Gryzinski,  one  of  which  is  Ref  16,  and 
Ochkur  and  Petrun’kin  (Ref  17).  Bauer  and  Bartkly  state 
that  their  results  agree  with  experimental  data,  where 
available,  within  a factor  of  3- 

Equations  20  and  21  are  two  solutions  they  find 
agreeable  with  experimental  results.  Equation  20  is  known 
as  the  Thompson  formulation  and  equation  21  as  the  Gryzinski 
formulation. 


/ H re4  \ 

Q1(E2>  ■(  \~)  % 


where 


V<e2-V 

(E2+l)(E2+l-Dn)I 


U <E0<U  , 

n-  2—  n+1 


Si  = 


Dn  <WU„> 
£ 


(E2+I)(E2+I-UnKE2+I-Un+1) 


, E2— Un+1 


(20a) 


(20b) 


21 


^2^E2^  ~ 


(21a) 


(Neire^)  g2  (E2,I,U)  6 ’ U, 


where  g9  = 


In 


E,+I-U 

2 n 

k 

/E2+I-Un\ 

T 

I 

3 

V V1  / 

/ /ul 

I o no  . 1 / _2 

\ 

1 I 

yd.  f d — 

/ 

(E2+i-V3  (E2+I)  . 

(21b) 


and 


3/2 


’ E2+I  ’ 

f un 

EQ+2I 

E2+I 

( — ± ) 
' 21+Ep-U 
2 n 


6 ' U = 


VUn 


U <E,<U 

, n—  2—  n+1 


Vl"Un  , 


E2-Un+1 


In  equations  20  and  21,  E2  is  the  energy  (eV)  of  the 
incident  electron,  I(eV)  is  the  ionization  potential  of  the 
state  from  which  the  atom  is  excited,  U (eV)  and  U , . (eV) 
are  the  bottom  and  top  energies  of  the  excited  level,  and 
Ng  is  the  equivalent  number  of  electrons. 

The  equivalent  number  of  electrons  (Ne)  is  an 
adjustable  parameter  that  can  take  on  several  values,  de- 
pending on  the  process  in  question.  For  electron  exchange 
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processes,  Bauer  and  Bartkly  recommend  Ng  equal  to  5. 


Bauer  and  Bartkly  also  identify  some  of  the  problem 
areas  of  equations  20  and  21 . The  function  in  equation  20 
fits  better  in  the  lower  energy  region,  but  drops  off  more 
rapidly  than  equation  21  in  the  high  energy  region.  Either 
equation  compares  reasonably  well  with  experiment,  and  the 
authors  feel  the  results  will  always  be  within  a factor  of 
4 of  experimental  data. 

Electron  Impact  Ionization 

The  electron  ionization  cross  section  formula  developed 
by  Gryzinski  is  used  throughout  the  Bauer  and  Eartkly 
article.  They  feel  the  Gryzinski  equation  agrees  best  with 
experimental  data.  The  formula  is  shown  in  equation  22. 

E 

QI  = Ne  gl  (22a) 

where 

1 V 1 3/2 

gl(X)  - 

(22b) 

ln(  2 . 72  + /FI)  , X = 

In  equation  22,  E2  is  the  energy  of  the  incident 
electron,  and  l(eV)  is  the  energy  necessary  for  ionization. 

Vrien  (Ref  18)  develops  a much  simpler  ionization 
cross  section  equation  that  compares  very  well  with  ex- 
perimental data.  The  results  are  shown  in  equation  23. 
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(23) 


740™p*  , 1.  1.  2J,  1 1_  ) 

yI  E+3.25I  V I E ' E ' j.2  E2  ; 

In  equation  23,  E is  the  incident  electron  energy, 

2 

I(eV)  is  the  energy  necessary  for  ionization  and  traQ  equals 
2 

0.88  X 10'  cm  . This  formula  works  best  when  the  energy 
for  ionization  is  less  than  3*5  ev. 

Equations  20  through  23  are  all  the  equations  necess- 
ary for  determining  pumping  rates  by  the  numerical  solu- 
tion of  the  electron  Boltzmann  equation  of  Chapter  1 . The 
next  section  calculates  the  cross  sections  for  Mercury  and 
compares  them  with  experimental  data,  when  available.  The 
cross  sections  eventually  used  in  the  numerical  calculations 
of  the  discharge  in  Chapter  4 are  identified  in  the  summary 
section  of  this  chapter. 

Results 

Rockwood  (Ref  5)  and  Cayless  (Ref  19)  show  the  cross 
1 3 

sections  for  the  6 Sq — ^*2  1 0 i;ransi‘kion • The  results 
are  illustrated  in  Figures  6 and  7,  respectively. 
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The  pumping  rate  constant  for  a specific  collisional 

process,  when  weighed  by  the  velocity,  is  proportional  to 

the  area  under  the  curve.  It  is  evident  that  the  6 Sq — 6^P2 

cross  section  is  several  times  greater  than  the  cross 

1 3 

section  for  the  6 Sq — 6Jp ^ q transition.  This  is  the  reason 
1 3 

that  only  the  6 Sq — 6JP2  ‘transi‘ti°n  is  tracked.  To  be 
completely  correct  the  other  transitions  should  be  included, 
along  with  the  reactions  among  the  6^p  . n level.  (See 

Appendix  A)  . 

The  Gryzinski  equation  of  the  previous  section  is 

1 3 

used  to  calculate  the  6 Sq — 6JP2  transition.  The  results 
are  shown  in  Figure  8.  Note  the  fall  off  at  high  energy 
which  may  cause  an  underestimation  in  the  total  pumping 
rate.  An  equivalent  number  of  electrons  of  20  is  needed  to 
bring  the  curve  up  to  the  same  order  of  magnitude  as  the 
curves  of  Rockwood  (Ref  6)  and  Cayless  (Ref  19).  This  value 
is  suspect,  but  no  other  method  is  available  for  a comparison. 
At  worst,  the  results  are  correct  to  an  order-of-magni tude . 
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Fig.  8.  Calculated  cross  section  from  Gryzinski  (G) 
and  Rockwood  (R). 
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The  6 Sq — 6 transition  of  Rock-wood  (Ref  5) 
and  Cayless  (Ref  19)  are  shown  in  Figures  6 and  9, 
respectively.  


20  30 

Electron  voltage  (V) 


Fig.  9.  Excitation  cross-section  for  for  Hg. 

A comparison  with  equations  19  and  20  of  the  previous 
section  require  an  equivalent  number  of  electrons  (Ng)  of 
20  to  bring  the  calculated  curve  within  an-order-of 
magnitude  of  Figure  9* 

The  6^-Sq — transition  cross  section  is  determined 

by  Yavorski  (Ref  12)  and  used  by  Cayless  (Ref  15)-  This 
cross  section  is  shown  in  Figure  10. 


jbVA 


Fig.  10. 
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Electron  voltage  (V) 


Cross-sections  for  excitation  of  7 from 
the  ground  state  and  6^P Q of  Hg,  after  Yav 
(Origin  transposed )( Ref  15)- 


Yavorski . 


3 3 

A collisional  cross  section  for  the  — 7 
transition  is  shown  in  Figure  11.  The  curve  is  a result 
of  the  work  of  Korotkov  and  Kazakov  - furnished  by 
Garscadden  (Ref  17)* 

Cross  sections  for  the  remaining  excitation  transi- 
tions are  not  found  in  the  literature.  The  Gryzinski 
solution  is  used  and  the  results  shown  in  Figures  12 
through  17. 

Any  pumping  rate  determined  from  these  collisional 
excitation  cross  sections  can  be  incorrect  by  an  order  of 
magnitude . Part  of  this  is  due  to  the  behavior  of  the  cross 
section  in  the  high  energy  region.  If  the  cross  sections 
for  transfers  within  the  triplet  or  singlet  region  of  the 
energy  level  diagram  behave  as  the  cross  sections  calculated 
by  Korotkov  and  Kazakov  (Ref  17)*  the  predicted  cross 
section  calculated  by  Gryzinski  is  too  low  in  the  higher 
energy  region.  Until  experimental  data  can  determine  these 
excited  cross  sections,  a full  understanding  of  the  micro- 
scopic processes  is  highly  speculative. 
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Rockwood  (Ref  5)  uses  the  ionization  cross  section 
shown  in  Figure  6.  This  result  is  very  similar  to  the  ioni- 
zation cross  section  used  by  Cayless  (Ref  19).  Figure  18 


shows  the  ionization  cross  sections  using  the  Gryzinski 
equation  and  the  Vrien  equation. 

Ionization  cross  sections  of  higher  metastable  levels 
of  Hg  are  not  available  in  the  literature.  Their  values 
are  calculated  using  the  Vrien  equation  and  shown  in  Figures 
19  through  22. 

The  last  section  of  this  chapter  is  a listing  of  the 
cross  sections  used  in  the  numerical  solution  of  the 
electron  Boltzmann  equation,  and  the  reference  for  the 
evaluation. 

Summary 

Cross  sections  not  included  in  the  following  list  can 
be  added  when  experimental  results  show  a need  for  their 
contributions . 

Table  II  lists  the  reactions  utilized  by  the  numerical 
HgCl  discharge  program  of  Chapter  4.  The  list  also  provides 
the  reference  for  determining  the  cross  section. 
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Fig.  19.  Ionization  of  Mercury  by  electron  impact. 
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Cross  Se 


Fig.  20.  Ionization  of  Mercury  by  electron  impact 
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Fig.  22.  Ionization  of  Mercury  by  electron  impact. 
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Table  II 

Mercury  Cross  Sections 


L 


Reaction  Reference 


6\— 

6\ 

5 

i\- 

6‘pl 

5 

6'Pi- 

6% 

16 

6lso~ 

73S1 

19 

6\- 

73S1 

16 

6\~ 

6\ 

16 

6\— 

73S1 

20 

6\- 

6\ 

16 

73S1  — 

6\ 

16 

63P2- 

6\ 

16 

&\— 

Hg+ 

5 

63P2— 

Hg+ 

18 

6lPl~ 

Hg+ 

18 

63S1  — 

Hg+ 

18 

Hg+ 

18 

The  following  two  sections  list  the  collisional  cross 
sections  for  the  element  Ar  and  the  dissociative  attachment 
cross  section  Cl2. 

Argon  (Ar) 

The  cross  sections  for  excitation  and  ionization  of 
mercury  are  well  documented.  These  values  are  widely  used 
and  considered  fairly  accurate.  The  listing,  and  Figures 
23  through  28  show  these  values. 
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Fig.  23.  Excitation  of  Argon  by  electron  impact. 
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Transition 

Energy  (eV) 

Reference 

Figure 

* 

At — Ar 

11.7 

21 

23 

** 

Ar — Ar 

13.2 

21 

24 

* ** 

Ar  — Ar 

1.5 

6 

25 

Ar — Ar+ 

15.759 

21 

26 

* + 

Ar  — Ar 

4.059 

18 

27 

+ 

Ar  — Ar 

2.559 

18 

28 

Chlorine  (Clg) 

The  dissociative  attachment  rate  for  Clg  is  determined 
by  treating  the  cross  section  as  a narrow  function  of  low 
energy.  The  value  of  this  function  is  9 X 10-1^cm2 
around  zero  energy  and  zero  everywhere  else.  This  tech- 
nique is  supplied  by  Hunter  (Ref  21). 

Conclusion 

The  results  accumulated  in  this  chapter  are  as  complete 
a list  as  can  be  found.  Several  of  the  cross  section  values 
have  been  determined  to  be  very  accurate,  while  others  vary 
up  to  a factor  of  9 in  comparison  with  the  same  cross 
section  determined  by  some  other  means. 

The  lists  are  not  meant  to  be  complete.  Their  in- 
tention is  to  give  an  order-of-magnitude  figure  for  de- 
termining the  temporal  evolution  of  the  HgCl  discharge. 

When  experimental  data  becomes  available,  they  should  be 
updated  with  this  new  information. 
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IV.  The  HgCIModel 


Overview 

In  this  chapter,  the  discharge  analysis  of  the  HgCl 
laser  is  continued  from  the  work  initiated  by  Summers 
(Ref  4).  The  chapter  is  divided  into  three  sections. 

First,  an  introduction  to  the  reactions,  second,  a simple 
analytic  approach,  and  third,  a discussion  of  results  from 
a complete  numerical  solution  of  the  HgCl  discharge. 

Introduction 

In  order  to  build  larger  scale  lasers,  one  must  de- 
velop a comprehensive  model  that  explicitly  accounts  for 
all  significant  processes.  Summers  (Ref  4)  recognizes  this 
need  and  identified  many  of  the  energy  flow  pathways.  He 
also  points  out  that  two  of  his  assumptions  would  need  some 
modification.  First,  the  assumption  that  the  secondary 
electrons  behave  as  a Maxwellian  distribution,  and  second, 
that  the  electron  beam  ionization  cross  section  for  producing 
secondary  electrons  is  not  a function  of  energy. 

The  first  assumption  is  eliminated  by  the  numerical 
solution  to  the  electron  Boltzmann  equation  identified 
in  Chapter  I.  The  second  assumption  is  adjusted  with  a 
technique  developed  by  Green  and  Sowada  (Ref  22). 

The  particular  species  which  control  the  processes 
in  the  HgCl  laser  are  electrons,  the  gound  and  metastable 
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levels  of  Hg,  the  ground  and  raetastable  levels  of  Ar,  and 
the  Cl2  molecule.  The  continuity  equations  of  the  dominant 
discharge  species  are  based  on  the  processes  listed  in 
Table  III.  The  table  is  neither  complete  nor  does  it  re- 
present a closed  system  of  reactions.  Equation  numbers 
identified  with  each  reaction,  and  typical  values  of  reac- 
tion rate  constants  determined  by  the  numerical  code,  are 
listed  in  Table  III. 


Table  III 

Processes  Included  in  Discharge  Model 


Reactions 

Rate/Cross 

Section 

Reference 

e-beam  ionization 

P + P 

e + Ar-*Ar  + e + e 

°1 

= 

7.6E 

! X 

10"l8cm2 

23 

(23) 

P + P 

e + Hg-*Hg  + e + e 

a2 

= 

8.1? 

’ X 

10~l8cm2 

22 

(24) 

electron  loss 

+ * 

e + Ar2  Ar  + Ar 

R1 

= 

3-9 

X 

1 0_8cm^/sec 

4 

(25) 

e + C12-^C1~  + Cl 

r2 

= 

1.0 

X 

10_^^cm^/sec 

Code 

(26) 

. -B-tt 

e + Hgg-^  Hg  + Hg 

R3 

= 

1.0 

X 

10-^cm^/sec 

11 

(27) 

metastable  formation 

* 

e + Ar  -»  Ar  + e 

% 

= 

1.0 

X 

10  ^2cm^/sec 

Code 

(28) 

* 

e + Ar  -*  Ar  + e 

R5 

= 

1.0 

X 

8 3 

10  crrr/sec 

Code 

(29) 

e + Ar  -*■  Ar  + e 

r6 

= 

1.0 

X 

10-^2cm^/sec 

Code 

(30) 

e + Ar  — - Ar  + e 

R7 

= 

1.0 

X 

8 3 

1 0 crrr/sec 

Code 

(3D 

-ft# 

e + Ar  -*  Ar  + e 

R8 

= 

1.0 

X 

10~^cm^/sec 

Code 

(32) 

-Mr 

e + Ar  -‘Ar  + e 

R9 

= 

1.0 

X 

6 3 

10'  cm;/sec 

Code 

(33) 
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Table  III  (Cont’d) 


Reactions 


Rate/Cross  Section 


Reference 


metastable  formation 


e 

+ 

Hg-»Hg  + e 

R10 

* 

e 

+ 

Hg  — * Hg  + e 

RU 

** 

e 

+ 

Hg  -*Hg  + e 

R12 

** 

e 

+ 

Hg  Hg  + e 

R13 

*** 

e 

+ 

Hg-*Hg  + e 

Rl4 

*** 

e 

+ 

Hg  — - Hg  + e 

R15 

e 

+ 

Hg^Hg  + e 

Rl6 

e 

+ 

Hg  — - Hg  + e 

R17 

* ** 

e 

+ 

Hg  Hg  + e 

Rl8 

**  * 

e 

+ 

Hg  — * Hg  + e 

R19 

e 

+ 

Hg  — * Hg  + e 

R20 

***  -* 

e 

+ 

Hg  -*  Hg  + e 

R21 

e 

+ 

Hg  Hg  + e 

R22 

****  * 

e 

+ 

Hg  — *■  Hg  + e 

R23 

e 

+ 

Hg  — * Hg  + e 

R24 

e 

+ 

Hg  — - Hg  + e 

r25 

e 

+ 

Hg  -*  Hg  + e 

R26 

e 

+ 

Hg  Hg  + e 

R27 

e 

+ 

Hg  — -Hg  + e 

R28 

****  -*** 

e 

+ 

Hg  — -Hg  + e 

R29 

ionization  formation 

e 

+ 

Ar-^Ar+  + 2e 

R30 

* + 

e 

+ 

Ar  — - Ar  + 2e 

r31 

1.0 

X 

10~^cm3/  sec 

Code 

(3*0 

1.0 

X 

10_®cm3/sec 

Code 

(35) 

1.0 

X 

10_^cm3/sec 

Code 

(36) 

1.0 

X 

8 o 

lO'cm-ysec 

Code 

(37) 

1.0 

X 

10  ^cmVsec 

Code 

(38) 

1.0 

X 

8 3 

10'°cmvsec 

Code 

(39) 

1.0 

X 

10  ^cm3/sec 

Code 

(40) 

1.0 

X 

10  cnr/sec 

Code 

(41) 

1.0 

X 

—8  3 

10  cnr/sec 

Code 

(42) 

1.0 

X 

10~^cm3/ sec 

Code 

(43) 

1.0 

X 

10  ^cm3/sec 

Code 

(44) 

1.0 

X 

-R  3 

10  cnr/sec 

Code 

(45) 

1.0 

X 

10-^cm3/sec 

Code 

(46) 

1.0 

X 

10~^cm3/sec 

Code 

(47) 

1.0 

X 

10-^cm3/sec 

Code 

(48) 

1.0 

X 

“8  3 

1 0 cnr/sec 

Code 

(49) 

1.0 

X 

—8  3 

10  cnr/sec 

Code 

(50) 

1.0 

X 

lO'^cm^/sec 

Code 

(51) 

1.0 

X 

n o 

10  'cnr/sec 

Code 

(52) 

1.0 

X 

lO'^cm3/ sec 

Code 

(53) 

1.0  X 10"15cm3/sec 

Code 

(54) 

1.0  X 10-Rcm3/sec 

Code 

(55) 
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Table  III  (Cont'd) 


Reactions  Rate/ Cross  Section  Reference 

ionization  formation 


**  + 

e + Ar  — *•  Ar  + 2e 

R32 

= 

1.0 

X 

lO'^cm-V  sec 

Code 

(56) 

e + Hg—  Hg+  + 2e 

R33 

= 

1.0 

X 

lO-^  cm^/sec 

Code 

(57) 

e + Hg*— Hg+  + 2e 

R34 

= 

1.0 

X 

—8  3 

10  cnr/sec 

Code 

(58) 

4. 

e + Hg  — Hg  + 2e 

R35 

= 

1.0 

X 

—8  3 

10  cnr/sec 

Code 

(59) 

***  + 

e + Hg  Hg  + 2e 

R36 

= 

1.0 

X 

lO^cm-Vsec 

Code 

(60) 

****  + 

e + Hg  — *>  Hg  + 2e 

R37 

= 

1.0 

X 

10~^cm^/sec 

Code 

(61) 

excimer  formation 

Hg+  + Cl -HgCl* 

R38 

= 

1.0 

X 

6 3 

10  cnr/sec 

4 

(62) 

Hg*  + Cl2— HgCl*  + 

Cl 

R39 

= 

3.8 

X 

1 sec 

11 

(63) 

Hg  + Cl2-  HgCl  + 

Cl 

R40 

= 

2.0 

X 

10  ^cmVsec 

11 

(64) 

***  * 
Hg  + Cl2—  HgCl 

+C1 

R4l 

= 

4.0 

X 

10  ^cmVsec 

11 

(65) 

****  * 
Hg  + Clg-*  HgCl 

+C1 

R42 

= 

4.0 

X 

1 0~-^cm-Vsec 

11 

(66) 

Hg2+  + Cl- HgCl*  + 

Hg 

r43 

= 

1.0 

X 

8 3 

10  cnr/sec 

11 

(67) 

Ar+  + Cl-  ArCl* 

r44 

= 

2.13  X 

[ 10-^cm-^/sec 

4 

(68) 

Ar*  + Clg— ArCl*  + 

Cl 

r45 

= 

3.16  X 

! 10_^®cm^/sec 

4 

(69) 

excimer  loss 

HgCl*—  HgCl  + hv 

R46 

a 

2.2 

X 

10^/sec 

11 

(70) 

* 

HgCl  + Cl2-»  Quench 

r47 

a 

4.6 

X 

10  ^cm^/sec 

11 

(71) 

* 

HgCl  + Ar  — Quench 

R48 

= 

9.0 

X 

10~*^cm-V  sec 

11 

(72) 

HgCl*  + Hg  — 2Hg  + Cl 

R49 

a 

1.0 

X 

10  ^cm-Vsec 

11 

(73) 

* 

ArCl  Ar  + Cl  + hv 

R50 

a 

5.0 

X 

lO^/sec 

4 

(74) 
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Table  III  (Cont'd) 


Reactions  Rate/Cross  Section  Reference 

3 body  reactions 


Ar+  + 2Ar-*Ar2+  + Ar 

R^  = 2.5  X 10~3^cm^/sec 

4 

(75) 

Cl  + Cl  + Cl2-*  2C12 

R^2  = I.165  X 10_32cm^/sec 

4 

(76) 

Cl  + Cl  + Ar-  Cl2  + 

Ar 

R53  = 1.957  X 10~32cm6/ sec 

4 

(77) 

Hg+  + Hg  + Ar  — Hg2  + 

Ar 

R^  = 1.0  X 10“31cm^/sec 

11 

(78) 

Hg  + Cl  + Ar  — HgCl  + 

Ar 

R^3  = 7-5  X 10  3^cm^/sec 

11 

(79) 

lower  level  loss 


2HgCl-Hg  + HgCl2  R56  = 3-0  X 10~10cm3/sec  11  (80) 

2HgCl-Hg2Cl2  R57  = 2.0  X 10~10cm3/sec  11  (8l) 

HgCl  + Cl2^HgCl2  + Cl  R^g  = 1.0  X 10-11cm3/sec  11  (82) 

Penning  Ionization 

Hg  + Ar+  + Ar  + e + Hg+  R^  = 4.37  X 10_12cm3/sec  Cal  (83) 
metastable  transfer 

Hg  + Ar-Hg  + Ar  R6q  = 1.0  X 10  i;cm7sec  11  (84) 
Hg  + Ar—Hg  + Ar  R6l  = 1.0  X lO"1  W/sec  11  (85) 

Rate  constants  identified  by  Tang  (Ref  11)  in  Table 
III  are  reasonably  accurate  though  uncertainties  of  rv*  50 % 
are  possible.  Reaction  rate  constants  calculated  by  the 
code  are  accurate  to  an  order  of  magnitude . 

In  Table  III,  the  reference  to  code  indicates  this 
reaction  rate  constant  is  dependent  upon  parameters  placed 
on  the  discharge.  Values  of  these  rates  are  calculated  by 


r 
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the  numerical  solution  of  Chapter  I and  are  not  listed  here. 
The  cross  sections  for  electron  beam  ionization,  reactions 
23  and  24,  and  the  calculation  of  Penning  ionization, 
reaction  83,  are  calculated  in  Appendix  B.  The  quenching 
products  of  reactions  ?1  and  72  are  considered  insignificant 
to  the  overall  gas  discharge. 

Table  III  is  as  complete  a table  of  reactions  as  can 
be  found.  Any  reaction  omitted  is  due  to  either  the  non- 
availability of  cross  section  information  or  the  reaction 
is  thought  to  be  insignificant  to  the  overall  process. 

In  the  next  section,  a simple  discharge  analysis  is 
developed.  The  purpose  of  the  section  is  to  illustrate  the 
physical  character  of  the  processes  taking  place.  For 
simplicity,  only  a few  of  the  reactions  of  Table  III  are 
used. 


Theory 


General  Reactions.  A rate  equation  can  be  generalized 


as  equation  86  where  n^  is  the  number  density  of  species  i. 

’ _ y / all  processes  \ _ y r all  processes  . r q&\ 
i ' creating  n^  ' ' annihilating  n.  ' ' ' 


annihilating  n^ 


A simple  example  of  some  of  these  processes  is  shown 
in  equation  87,  where  i represents  the  upper  laser  level. 


/ electron  impact  s / chemical  •, 

pumping  ' ' pumping  ' 


, 


/ stimulated  _ , spontaneous 
' emission  ^ ~ 1 emission 
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In  equation  87 » the  term  electron  impact  processes 
are  those  processes  which  create  or  destroy  metastable 
levels  by  collisions  of  the  first  or  second  kind  (Ref  24). 
The  term  stimulated  emission  deals  with  optical  excitation 
and  deexcitation,  v/hile  the  term  spontaneous  emission  is  the 
process  of  radiative  relaxation.  A full  understanding  of 
the  chemical  pumping  term  is  beyond  the  scope  of  this  paper. 

The  next  section  is  the  development  of  the  discharge 
analysis.  A complete  analysis  is  virtually  impossible, 
therefore,  reactions  that  characterize  the  overall  physics 
of  the  process  are  used. 

Model 

The  rate  equations  88  through  90  are  assumed  to  be 
the  controlling  processes  in  the  HgCl  discharge  evolution. 

Electrons 


n 


= S,  - R2nencl^ 


R33nenHg 


R34nenHg* 


(88) 


* 


^g*  = ^oV^fe  ~ RllnenHg*  + R19nenHg** 


(89) 


'34nenHg*  " R39nHg*nCl?  " Rl8nenHg* 


SI 2 

nCl2  = " R2nenCl2  " R39nHg*nCl2 


(90) 
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r 

In  equation  88,  3^  is  a source  term  of  secondary 
electrons  generated  from  the  electron  beam.  A development 
of  these  source  terms  is  explained  in  Appendix  B.  Each 
term  in  equations  89  through  90  represents  a different 
reaction  of  Table  III.  The  reaction  process  for  each  of 
these  terms  is  identified  with  the  appropriate  rate  symbol. 

Equations  89  through  90  can  be  solved  to  determine 
both  the  temporal  evolution  of  the  discharge,  and  the  con- 
ditions for  the  existence  of  steady-state  solutions.  Once 
the  conditions  are  determined,  a study  of  the  stability  can 
be  made. 

The  rate  equation  for  the  electron  number  density, 
equation  88,  is  determined  by  a balance  between  ionization 
rates  and  electron  attachment  rates.  Steady  state  opera- 
tion is  never  quite  achieved  because  the  electron  number 
density  increases  as  the  rate  of  dissociative  attachment, 
equation  90,  decreases. 

Summers  (Ref  4)  identified  two  important  points  that 
aid  in  the  analysis.  First,  the  depletion  of  Cl2  is  slower 
than  the  rate  of  ionization  of  Hg  by  at  least  an  order-of- 
magnitude.  Second,  the  ground  state  neutral  species  vary 
on  time  scales  of  microseconds  while  ionized  species  vary 
on  the  order  of  nanoseconds.  Because  of  these  features, 
equilibrium  can  be  reached  o..  short  time  scales  even  though 
long  term  equilibrium  is  impossible. 

The  equilibrium  of  metastable  species,  equation  89, 
can  be  controlled  by  either  a dominant  eximer  formation  or 
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by  a dominant  ionization  of  metastable  levels.  It  is  of 
interest  to  examine  both  of  these  conditions. 

When  excimer  formation  is  dominant,  equation  89 
becomes  equation  91  • 


loVhte 

v r\j  ■ ■ ■ -w. 

Hs*  R39nci2 


(91) 


The  electron  rate  equation  now  becomes  equation  92. 


ne  [ R33"Hg  - R2nCl2  ^ + "e2  C ]<92) 


n = S.  + 
e 1 


In  equilibrium,  ng  = 0,  the  electron  number  density 
is  determined  by  solving  the  quadratic.  The  time  dependent 
solution  is  obtained  by  directing  integration. 

Fortunately,  Hunter  (Ref  21)  has  already  evaluated  a 
similar  function  with  the  KrF2  laser.  He  shows  that  the 
eximer  formation  equilibrium  is  a condition  corresponding 
to  an  electron  beam  dominated  discharge. 

If  the  electron  density  becomes  large  enough  for 
metastable  loss  rates  to  be  dominated  by  ionization, 
equation  89  becomes  equation  93* 


n * „ !!lOne_Hg.  = 

Hg  R34ne  R34 


(93) 


The  electron  density  rate  equation  is  now  approximated 
by  equation  94. 


J 
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ne  * S1  + *"«: 


In  equation  95,  y is  the  electron  attachment  rate 
minus  the  cumulative  ionization  rate. 

The  time  dependent  solution  is  found  by  integration, 
and  the  equilibrium  condition  is  simply 


(96) 


Hunter  (Ref  21)  also  shows  that  this  equilibrium 
corresponds  to  a field-dominated  discharge.  His  results 
are  shown  in  Figure  29* 

The  temporal  evolution  of  the  discharge  varies 
according  to  the  concentration  of  Cl2  and  Hg.  If  the  eximer 
formation  equilibrium  is  obtained,  the  depletion  of  Cl2  will 
cause  the  discharge  to  transition  to  ionization  equilibrium. 
This  is  caused  by  the  continual  depletion  of  Cl2  by  disso- 
ciative attachment.  Once  the  electrons  reach  a certain 
density,  the  depletion  of  Cl2  rapidly  increases  by  the  in- 
creasing secondary  electron  density  and  breakdown  is  forth- 
coming. The  production  of  electrons  in  the  ionization 

I 


6o 


Number  Density  (cm 


0 100  200  300  400 

time  (nanoseconds) 

Fig.  29  Temporal  evolution  of  KrF  discharge (Ref  21). 


equilibrium  region  is  approximately  equal  to  the  number 
density  of  ionized  mercury  because  of  its  low  ionization 
energy. 

Intensity 

Another  interesting  feature  to  evaluate  is  the  in- 
tensity rate  of  the  laser  discharge.  The  stimulated 
emission  .term  of  equation  87  can  be  rewritten  as  equation 
97  (Ref  26). 


An 


m < Nu  - ij  ) 


(9?) 


In  equation  97 » the  subscripts  u and  z represent  the 
upper  and  lower  laser  levels,  respectively.  The  N re- 
presents the  number  density  of  the  appropriate  level,  and 

-2  -1 

n is  the  photon  flux  in  cm  sec  . The  symbol  g stands 
for  the  degeneracy  of  that  state,  while  the  symbol  a is  the 
stimulated  cross  section  shown  in  equation  98  (Ref  26). 


a 


8tt2tc6X 


(98) 


where  X is  the  wavelength  in  cm 

. -1 
T is  the  Einstein  A coefficient 

c equals  the  speed  of  light  in  cm/sec 

6X  is  the  half  width  of  half  max 


Assuming  the  degeneracies  of  the  upper  and  lower  level 
to  be  equal,  the  gain  of  the  medium  is  simply  equation  99. 
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g = o [ Nu  - N,  ] 


(99) 


r 


In  general,  the  change  in 
medium  is  shown  in  equation  100 


intensity  of  a laser 
(Ref  25). 


I = cl  ( g - gth  ) 


(100) 


where  g^.^  represents  the  laser  threshold  gain.  Hunter 
(Ref  21)  derives  a more  complete  intensity  equation,  which 
include,  a term  due  to  the  spontaneous  radiation.  This  is 
shown  in  equation  101. 


d ( g - gth  ) 


+ 


chvN^dfi 

T spontaneous 


(101) 


In  equation  101,  hv  is  the  energy  of  the  photon,  dn 
the  solid  angle  based  on  the  length  of  the  cavity,  and  t 
the  inverse  Einstein  A coefficient. 

The  next  section  discusses  the  results  generated 
from  the  numerical  solution  of  the  HgCl  laser. 

Discussion 

It  is  interesting  to  note  some  of  the  results  of  the 
numerical  solution.  The  time  to  breakdown  is  very  sensitive 
to  the  concentration  of  Hg.  With  a very  small  field,  l 
Townsend  or  less,  a beam  of  1 amp/ cm  caused  breakdown 
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inside  of  100  nanoseconds,  with  a ratio  of  Ar/Hg/Clg  as 
95%/95?°/^  • &%/0 . 2%.  Reducing  the  concentrations  to  match 
those  of  Tang  (Ref  1),  breakdown  occurred  in  a much  longer 
time  period.  The  sensitivity  of  breakdown  to  Clg  is  ex- 
pected since  this  is  the  only  "sink"  term  for  electrons. 

It  is  also  interesting  to  note  that  the  addition  or 
removal  of  the  "stored  energy"  levels  of  the  mercury  atom, 
did  not  significantly  effect  the  temporal  evolution  of  the 
discharge . 


V.  Results  and  Recommendations 


Results 

The  purpose  of  this  paper  is  to  determine  the  temporal 
evolution  of  the  HgCl  discharge.  It  has  reached  that  objec- 
tive. Before  presenting  the  results,  a few  words  about  the 
concluding  remarks  are  apropos. 

It  is  always  easiest  to  arrive  at  a conclusion  when 
there  are  experimental  data  to  back  the  results.  In  the 
case  of  the  HgCl  laser,  there  aren't  many  experimental  data, 
but  there  is  enough  to  give  credibility  to  any  matching 
numerical  solution. 

The  results  of  this  paper  are  determined  more  from 
basic  physics,  and  knowing  a "ball  park"  figure,  than  from 
any  great  esoteric  derivation.  They  are  by  no  means  an 
attempt  to  discredit  any  published  literature.  There  are 
three  results  that  point  out  areas  that  need  more  research 
in  order  to  better  model  the  HgCl  laser. 

First,  a complete  numerical  solution  of  the  electron 
Boltzmann  equation  was  determined.  This  numerical  solution 
was  checked  against  two  sources  of  information  and  agreed 
remarkably  well. 

When  fractional  energy  deposition  was  compared  to 
Rockwood  (Ref  5)>  a definite  difference  was  noted  in  the 
lower  E/N  values  (Figure  2).  Rockwood  observed  a bend  in 
this  region,  while  the  results  of  the  numerical  solution 


used  in  this  paper  did  not  match  those  results.  This  was 
first  thought  to  be  the  effect  of  a coarse  bin  width  of 
1 eV.  After  including  elastic  collisions,  it  was  determined 
that  the  bend  was  the  result  of  elastic  processes  dominating 
the  discharge.  The  original  assumption  of  considering 
elastic  collision  processes  to  be  insignificant  had  to  be 
abandoned  if  results  were  to  be  reliable  in  all  ranges 
of  E/N. 

Second,  the  cross  sections  for  transitions  among 
metastable  levels  of  Hg  are  calculated  with  tested  known 
theory.  There  appears  to  be  a small  difference  between 
calculated  and  published  information  on  cross  sections,  how- 
ever, a more  significant  point  is  the  sensitivity  of  the 
reaction  rate  constants  of  these  cross  sections.  Varying 
a cross  section  by  an  order  of  magnitude  varies  the  reaction 
rate  constant  by  at  least  an  order  of  magnitude.  It  is  also 
found  that  the  overall  discharge  characteristic  did  not 
change  when  including  the  higher  stored  energy  levels.  The 
most  significant  contribution  to  the  discharge  is  the  ion- 
ization  of  the  6^Pg  level  of  Mercury.  As  soon  as  this  level 
became  significantly  populated  (~10  -^crn  J),  breakdown  oc- 
curred in  a matter  of  nanoseconds.  This  is  due  mostly  to 
the  ionization  process  of  the  level. 

Third,  the  temporal  evolution  of  the  HgCl  laser  yields 
results  that  are  somewhat  interesting.  This  third,  and  most 
significant  feature,  calls  for  some  speculation  on  the  in- 
formation supplied  to  the  numerical  HgCl  laser  program. 
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On  the  initial  runs  of  the  HgCl  laser  program, 
secondary  electron  number  densities  were  building  up  in  a 
remarkably  short  time.  With  no  E field  on  the  gas  and  an 
electron  beam  of  1 amp/cm  , breakdown  occurred  in  less  than 
100  nanoseconds.  A figure  of  400  to  500  nanoseconds  is  more 
reasonable . 

The  rapid  buildup  of  secondary  electrons  was  first 
thought  to  be  the  result  of  ignoring  collisions  of  the 
second  kind  - it  wasn't.  Perhaps  the  idea  of  stored  energy 
in  the  higher  metastable  levels  of  mercury  over-estimated 
the  excited  number  densities.  This  adjustment  did  not  make 
an  appreciable  change.  Finally,  a reduction  in  the  electron 
exchange  cross  sections  gave  the  results  more  physical 
meaning.  Gas  breakdown  did  not  occur  up  to  times  of  250 
nanoseconds  or  greater.  The  gas  reached  a transient  stable 
condition.  By  increasing  either  the  electron  beam  or  the  E 
field,  breakdown  occurred  in  shorter  time  periods. 

It  is  also  evident  that  the  concentration  of  mercury 
strongly  controls  the  electron  density.  A factor  of  5 re- 
duction in  the  concentration  increased  the  time  to  break- 
down significantly. 

The  reactions  of  Table  III  are  the  ones  that  determine 
the  twenty-one  differential  equations  to  be  solved.  The 
numerical  solution  of  these  differential  equations  yield  the 
temporal  evolution  of  the  HgCl  discharge.  Although  the 
stiffness  of  the  differential  equations  prohibited  the  in- 
tegrator from  solving  the  equations  beyond  200  nanoseconds, 
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a general  trend  can  be  seen  as  compared  to  the  KrF  discharge. 
This  is  shown  in  Figure  30. 

Small  variations  in  the  initial  concentration  of 
CI2  also  effect  the  electron  build-up  significantly.  The 
concentrations  were  finally  adjusted  to  match  the  values 
used  by  Hunter  (Ref  1).  In  his  experiment  he  also  noted 
the  sensitivity  of  the  discharge  to  the  concentration  of  Cl2. 
It  is  the  above  results  that  lead  to  the  following  recom- 
mendations . 

Recommendations 

In  summary,  work  needs  to  be  continued  in  three  major 
areas.  First,  a closer  examination  of  the  excitation  cross 
sections  for  Hg  should  be  attempted.  This  can  be  accom- 
plished by  either  redoing  the  work  of  Rockwood  (Ref  5)  or 
by  attempting  a quantum  mechanical  calculation  similar  to 
that  of  Yavorsky  (Ref  12).  Second,  a determination  of 
several  other  reactions,  should  be  evaluated.  This  evalua- 
tion should  include  an  examination  of  the  lower  levels 
of  mercury.  Nikitin  (Ref  27)  suggests  that  the  -^P0  level 
may  become  heavily  populated  from  the  higher  levels  of 
mercury.  The  examination  should  also  include  a means  of 
removing  the  lower  levels  of  HgCl.  Burnham  (Ref  28)  shows 
a marked  increase  in  efficiency  with  the  addition  of  the 
diatomic  molecule  N2.  Third,  experimental  comparison  must 
be  continued.  This  part  can  best  be  accomplished  by  finding 
a better  integrator  for  the  stiff  differential  equations. 
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Number  Density  (cm 


1 


Conditions : 

E/N  ~ 0.0 

Ar  Concentration  = 
Hg  Concentration  = 
Clp  Concentration  = 
j = 1 amp/cm 

N^t  = 3 X 10l9cm 
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Fig.  30  Temporal  Evolution  of  HgCl  Discharge. 
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The  HgCl  laser  is  in  the  infancy  of  development  and 
the  pathwork  for  understanding  a large-scale,  viable  laser 


is  laid  out.  The  gain  of  insight  toward  the  microscopic 
processes  will  be  invaluable  in  any  e-beam  controlled  dis- 
charge laser. 


i 
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APPENDIX  A 


Transfers  among  the  6-^P2  ^ o s^a^es  3X6  no^  widely 

O O 

known.  Yavorsky  (Ref  12)  calculates  the  — 6^p2  tran- 

sition with  the  result  shown  in  Figure  31. 


Fig.  31.  Cross-sections  for  transitions  of  the 
first  and  second  kind  between  6-^P.  and 
6^P2,  after  Yavorski . 

Unfortunately  no  data  on  the  other  transitions  exists. 
Cayless  (Ref  19)  uses  a simple  Boltzmann  population  argument 
to  show  that  the  cross  sections  for  excitation  from  o e 
state  to  another  does  not  differ  greatly  from  one  another. 
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APPENDIX  B 


Penning  Ionization 

— IP  P 

The  reaction  83  of  4.37  X 10  cnr/sec  is  determined 
by  the  equation 


R 


<c rv> 


m o(€)61//2 f(€)  d€ 

OO 

/ f(€)  d€ 


(1) 


2 

where  € is  energy  (eV),  a is  cross  section  (cm  ) and 

f(€)  is  the  number  of  atoms  per  cnr  per  eV.  Assuming  the 

* 

distribution  of  heavy  species  (Hg  + Ar  ) are  in  a Maxwellian 
Distribution  equation  (1)  becomes  equation  (2) 


R = av  = a 


8kT 

irMr 


1/2 

) . Mp  = 


M1  M2 
+ M2 


T = 300°K  K = Boltzmann  constant 
M1  = 200.59  AMU  (Hg) 

M2  = 39.95  AMU  (Ar) 


(2) 


Using  the  cross  section  determined  by  Phelps  (Ref  30) 
_iA  2 

of  1 X 10  cm  for  Penning  ionization  results  in 


R = 4.37  X 10'12cm3/sec 
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e-beam  ionization 


The  approach  used  in  calculating  the  reaction  rate 
constant  is  the  one  utilized  by  Green  and  Sawada  (Ref  22) 
and  by  Bass,  Berg,  and  Green  (Ref  23). 


S(E,T)  = A(E) 


-HE)2 

[T-T0(E)]2  + T (E)2 


(1) 


where 


A(E)  = 10"l6cra2  | In  ( | ) , 


r(E)  = r 


s E + fb  ' T0(E)  ' Ts  ~ E + Tv 


and  the  units  of  S(E,T)  are  cm2/eV. 


o = S S(E,T)  dT  = total  ionization  (2) 

cross  section 

a = A(E)  (E)  [ tan"1  ( ^P---*0  ) + tan"1  (—)'}, 
where  \ ~ \ (E-I) 

The  parameters  used  for  Hg  are  those  of  Smith-Liska 
in  Ref  23,  while  those  for  Ar  are  the  ones  listed  in  Ref  22. 
An  Energy  (€)  of  300KeV  is  used  for  the  electron  energy.  The 
units  of  S(E,T)  and  Tffl  are  cm2/eV  and  eV  respectively. 
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APPENDIX  C 


Numerical  Code 

The  numerical  solution  of  the  Boltzmann  time-dependent 
equation  is  the  joint  effort  of  several  people.  By  project- 
ing the  equation  on  an  energy  axis  of  K tins  with  width  a?, 
the  result  is  a set  of  K coupled  differential  equations. 

The  form  of  these  equations  can  be  put  in  the  following 
notation  (Ref  5)* 


[Jj(k)  - Jf(k)]  [J+X(k)  - Jel(k)] 
A€  A€ 


+ s j Ns^Rsjk+m  .+Rsjk-m  .nk-m  .Ns//Rs+Rsk+m„.nk+m„. 

9 d S J SJ  SJ 


SI  SI 


with 


+ m Rsmnm“^Rs jk+Rs jk+Rsk^nk^ 


+ n , 2Ne2  , E ,Z  €k  , nk+l+nk  nk+l  nk  N 
JfU)  ~ ~3m~  ( N } v+7N(  ) * 


V"  4€k  A€ 


= "4  C ( f ~ ) -kT<  ( 


) 1 


A€ 


and 


vk  2€k  !/2 

TT“  <-“> 


s s ^ k; 


-+  0 zejvfe  qsas(e+) 

= 2mN  — E — k_ 
k m £ m 

s s 
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where 


+ — + 

€k  = kA€’  €k  = €k-l  ’ 

* 

J"(k)  = J+(k-l),  mgj  = €sjAe 

The  complete  code  is  maintained  ty  Captain  (Major 
Selectee)  Hunter  of  the  Physics  Department  at  AFIT. 
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